Cosmic relic abundance and f(R) gravity 
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The cosmological consequences of f(R) gravity are reviewed in the framework of recent data 
obtained by PAMELA (Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics) 
experiment. This collaboration has reported an excess of positron events that cannot be explained 
by conventional cosmology and particle physics, and are usually ascribed to the dark matter presence 
(in particular, weak interacting massive particles). The dark matter interpretation of PAMELA data 
has motivated the study of alternative cosmological models (with respect to the standard cosmology) 
owing to the fact that they predict an enhancement of the Hubble expansion rate, giving rise, in 
such a way, to thermal relics with a larger relic abundance. Our analysis shows that f(R) cosmology 
allows to explain the PAMELA puzzle for dark matter relic particles with masses of the order or 
lesser than 10 2 GeV in the regime p c < p m where p c is the curvature density and p m the radiation 
density. For the model f(R) — R + aR n , it then follows that n ~ 1 and small corrections with 
respect to General Relativity could lead indeed to address the experimental results. However other 
interesting cosmological models can be considered during the pre-BBN epoch as soon as the BBN 
constraints are relaxed. In such a case, the PAMELA data can be fitted for a larger class of f(R)- 
models. 

PACS numbers: 04.50.-h, 98.80.-k, 98.80.Es 



(N 
> 

o 

(N 

o 

(N 



X 



I. INTRODUCTION 

The discovery of the accelerated expansion of the Uni- 
verse [1] has motivated the developments of several alter- 
native theories of gravity. These models are built up, typ- 
ically, either in the framework of the conventional Gen- 
eral Relativity (GR) or in the framework of its possible 
generalizations or modifications. In the last years, among 
the different approaches proposed to generalize Einstein's 
GR, the so called f(R) gravity has received a growing 
attention. The reason relies on the fact that such an ap- 
proach allows to address the problem of the observed ac- 
celerating phase of the Universe, without invoking exotic 
sources of dark energy but just extending the standard 
paradigm of GR. The gravitational Lagrangian for these 
theories is a generic function of the Ricci scalar curva- 
ture R (not necessarily linear as in the Hilbert-Einstein 
action), and the corresponding action with the inclusion 
of standard matter reads 



S '= J dAx V-9f(R) + S m [g^,ipm] 



(LI) 



where k 2 = 8irG. It is a difficult task to deal with 
higher order terms in the scalar curvature, thus the forms 
of f(R) that have been most studied in literature are 
analytic functions of the forms f(R) ~ R + aR n or 
f(R) ~ R n with n positive or negative. Also broken 
power laws have been extensively considered [H-IH- 

Current cosmological observations indicate that our 
Universe is dominated by dark matter, responsible of 
galactic and extragalactic dynamics as well as of struc- 
ture formations, and dark energy, responsible of the ac- 



celerated expansion of the Universe. In other words, 
these unknown ingredients mean that we need a source 
to cluster large scale structure and a source to speed up 
the Hubble fluid. The ratio between the cold dark matter 
and the dark energy with the critical density satisfies the 
bounds [§] 

0.092 < QcDMh 2 < 0.124 , 0.30 < n DE h 2 < 0.46 , 

(1.2) 

where h = lOOKm s 1 Mpc 1 . Although many models 
have been proposed, the nature of the two components is 
still unknown. Favorite candidates for non-baryonic cold 
dark matter are the so called WIMPs (weakly interacting 
massive particles). However, such particles have not been 
detected yet. 

The aim of this paper is to explore the implications 
of the modified cosmology provided by f(R) gravity to 
the thermal relic abundance. Alternative cosmologies in- 
deed predict modified thermal histories for relic particles. 
These modifications, in principle, occur during the pre 
big bang nucleosynthesis (BBN) epoch, a period of the 
Universe evolution not directly constrained by cosmolog- 
ical observations. To account for the enhancement of the 
expansion rates that cosmological models provide with 
respect to GR, it is usual to write [l(| 



H(T) = A{T)H GR {T) 



(1.3) 



where H and Hgr are the expansion rates of the alter- 
native cosmological model and of GR, respectively. The 
function A(T) is the enhancement factor. Its temper- 
ature dependence is such that A(T) > 1 at large tem- 
perature and A(T) —> 1 before BBN set up, i.e. Eq. 



2 



(|I.3|) holds at early times, while at later time H = Hqr- 
This last requirement is necessary due to the successful 
prediction of BBN on the abundance of primordial light 
elements. Different cosmological scenarios have been 
proposed in literature to constrain the function A(T) 

pi 111 EMI- 



When the expansion rate of the Universe is enhanced 
(as compared to that one derived in the framework of 
GR), thermal relics decouple with larger relic abundance. 
The change in the Hubble rate may have therefore its 
imprint on the relic abundance of dark matter, such as 
WIMPs, axions, heavy neutrinos. This kind of studies is 
motivated by recent astrophysical results which involve 
cosmic ray electrons and positrons [191422} . antiprotons 
[23} , and 7-rays 0, [H| ■ Particular attention is devoted 
to the rising behavior of the positron fraction observed 
in PAMELA experiment Besides the astrophysical 
interpretation of this phenomenon (26|, it is under in- 
vestigation the possibility that the raising of A(T) could 
be due to dark matter annihilation, dominantly occurred 
into leptons [I3,[28j]. In this last case, it is required a large 
value of (a ann v). More specifically, PAMELA and ATIC 
(Advanced Thin Ionization Calorimeter) data require a 
cross section of the order or larger than (a an nv) ~ 10 26 
cm 3 sec -1 : in this case thermal relics would have the ob- 
served dark matter density [29]. In this paper, we want 
to investigate if such data and constraints can be framed 
in the context of f(R) gravity. 



The paper is organized as follows. In Sec. II we derive 
the f(R) gravity field equations and discuss the energy- 
momentum tensor related to the higher-order curvature 
terms. In particular, we emphasize the possibility to de- 
scribe the curvature induced terms as an effective per- 
fect fluid. In Sec. Ill we write down the expressions for 
the energy density and pressure induced by higher order 
curvature-terms for a power law function of the scale fac- 
tor. We then consider examples of f(R) models in which 
the ratio of the curvature energy density and curvature 
pressure is constant. In Sec. IV, we discuss the cosmic 
enhancement in view of the abundance of thermal relics 
recalling the main results reported in [Io[. Sec. V is de- 
voted to the study of the amplification factor given by 
(|I.3p in the framework of f(R) models. The pre-BBN 
epoch is considered in Sec. VI. In such a case, f(R)- 
models are less "fine-tuned" and several of them, in prin- 
ciple, could fit the PAMELA data. In Sec. VII we apply 
the conformal transformations to recast the f(R) gravity 
in terms of GR plus a scalar field. The latter gives rise 
to an effective potential related to the from of f(R) and 
its derivative. Conclusions are drawn in Sec. VIII. In 
Appendix we report some useful formulas adopted in the 
main text. 



II. FIELD EQUATIONS AND 
ENERGY-MOMENTUM TENSOR IN f(R) 
GRAVITY 

The variation of the action (jl.ip with respect to the 
metric yields the field equations 

f'R,* - { 9»u - V M V„/' + g^Uf = k 2 T™ , (II.l) 

where the prime indicates the derivative with respect to 
R and T™ is the energy-momentum tensor of matter. 
Here m stands for radiation (rad) or standard matter 
(mat). The trace reads 

3D/' + f'R - 2f = n 2 T m . (11.2) 

For our aim, it turns out convenient to rewrite (|II.1[) in 
the form 

g^ = k 2 (t;i + t^), (11.3) 

where G^ v is the Einstein tensor 

G[iV R^iV ^ 9jlvR ! 

and T^ v is the curvature energy-momentum tensor in- 
duced by higher order terms in the curvature invariants 
[l2T |. It is defined as 

k 2 T* u = (1 - f')R^ + \(f- R)9n» + V^uf - g^nf . 

(II.4) 

It gives rise to an effective description of the source term 
of Einstein field equations. A comment is in order. In 
the right-hand side of Eq. (|II.3p two effective fluids appear 
indeed: a curvature fluid and a standard matter fluid. 
This representation allows to treat fourth order gravity 
as standard Einstein gravity in presence of two effective 
sources [30(. This means that such fluids can admit fea- 
tures that could be unphysical for standard matter. Con- 
sequently all the thermodynamical quantities associated 
with curvature should be considered effective and not 
bounded by the standard constraints related to matter 
fields. Moreover, this description does not compromise 
any of the thermodynamical features of standard matter 
since Bianchi's identities are separately fulfilled for both 
fluids as we will show below. 

By using the properties [V 7 ,V^]V Q = — R a p/3j V p , 
where V a is a generic vector, and [V^V^]/' = 0, it 
is straightforward to show that 

V ai (k 2 T c ^) = 0. (II.5) 

As a consequence, the divergences of both sides of Eq. 
pi.3p vanishes, provided 

V M T m ^ = 0. (II.6) 

The latter equation implies that for f(R) gravity one can 
separately require that the energy-momentum tensor of 
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matter is preserved, independently as the gravitational 
background evolves. 

The aspect that arises from the above analysis is that 
one can define the total energy-momentum tensor as 

rp^Ll/ rjmi fJ,V _|_ rpC ^lV 

and then V M T Ml/ = 0. Moreover, it is possible to re- 
interpret the energy-momentum tensor induced by the 
curvature as a perfect fluid (we are assuming a homoge- 
neous and isotropic Universe) [30T ] 

t° " = (p c + p c Wu» - P c s; = ( P c , - P c , - P c , - P c ) , 

(II.7) 

where p c and p c are the energy density and pressure in- 
duced by the curvature terms; u^ 1 = (1,0) is the four- 
velocity of the effective fluid. For matter, treated as a 
perfect fluid, one has 

77" = ( P m +p m )u fi u v -p m s; = (fr,-p m ,-p m ,- P m ). 

(11.8) 

For a spatially flat FRW's metric 



ds 2 = dt 2 - a 2 (t)[dx 2 + dy 2 + dz 2 } , 



(II.9) 



one can write the explicit expressions for p c = k 2 T c ° and 
-p c 5\ = k 2 T\ j as 



Moreover, the equation of continuity for matter gives 



p" 1 = p^a~ 3(1+t ' 



(11.16) 

where w = p m /p m is the adiabatic index. We have de- 
liberately indicated the " adiabatic" indices a and w with 
different symbols in order to point out the different con- 
tributions to the cosmic dynamics of the two fluids. 



In what follows we will analyze the regime p 



c < p m^ 



i.e. the effective curvature is sub-leading with respect 
to matter. As we shall see, the alternative cosmology 
provided by f(R) gravity allows to explain PAMELA'S 
results for tiny deviations from GR. While the approach 
to matter fluid is standard, the curvature fluid requires 
a detailed discussion as discussed in the next section. 



III. POWER LAW EXPANSION OF THE 
COSMOLOGICAL BACKGROUND 

As well known, the high nonlinearity of f(R) gravity 
makes quite difficult to get exact solutions of field equa- 
tions 4]. Here we are interested in finding cosmological 
scenarios where a is constant in order to compare it with 
the proper matter adiabatic index w. 

A possibility to get a constant a is provided by a power 
law evolution of the scale factor 



P c = (i - f)R° + -(f-R)- mf , (ii.io) 



p c = -(l - f')R\ --(f-R) + f + 2Hf , (11.11) 



a(t) = a t p . 



(III.1) 



where 



f' = f"R, f' = f"'R 2 + f"R, 



and H — a /a is the expansion rate of the Universe. In 
the above formula, there is not sum over the indices i in 
R\ (the latter represent the diagonal components of the 
Ricci tensor), while the dot stands for the derivative with 
respect to the cosmic time t. 

From (|II.5|) one immediately gets 

p c + 3H(p c +p c ) =0. (11.12) 

In particular, setting 

p c = ap c , (11.13) 

the solution of (|ILT2|) assumes the standard form 

p c = p c e- 3 -f H ( 1+ ^ dt . (11.14) 

where Pq is a constant. 

In the case in which a is independent of cosmic time, 
as provided by some models of f(R) gravity, the previous 
equation reads 



Such a power law behavior holds in the regime in which 
we are interested, that is in presence of both curvature 
and matter fluids. 

Eqs. (|II.10|) and (|II.11|) (see Appendix for details) as- 
sume the form 

» c = -wh) R+ ^~ wV fl - 3H/ '' (m - 2) 

y 6(2/3-1) 2 6(2/3- FT J 3 ' 

(III.3) 

Notice that for f(R) = R, both p° and p c vanish, as ex- 
pected. Therefore corrections to the Einstein-Hilbert ac- 
tion provide a non-trivial expression for a. In Appendix, 
we present some details aimed to show that (|II.L2[) is 
explicitly fulfilled. 

It is worthwhile to write down the explicit expression 
of o as 



p c 1 Af 



p c 3D 

N = (3(3 - 2)R - 3(2/3 - 1)/ + (3/3 - l)f'R + 



2(2/3-3) , 2 4 3 
— f R ~~p f R ' 



(III.4) 



P C =pC a -^ 



(11.15) 



V = -/3R + (2/3 - 1)/ + (1 - f3)f'R - 2f"R 2 
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TABLE I: In the Table, expressions of a for different f(R) 
models are reported. The quantities c(/3, n) and d(/3, n) are 
constants which depend only on /3 and n (see Eqs. (IIII.5|l and 
(1HI6)) '). These values are independent of the adiabatic index 
w of matter due to the validity of Bianchi identities for matter 
and curvature fluids, which allows to write two independent 
equations of state. 





f(R) 


a 


1. 


R + aR 2 


4-3/3 
3/3 


2. 


R + aR" 


3d(/3,n) 


3. 


R + aR 2 + -yR s 


1 a(4 - 3/3) + 7(4 - p)R 
/3 3a + 7 (/3 + 10)R 


4. 


yR n 


13/3 - 2 + -yc(/3,n)R n - 1 
3 -/3 + -yd(/3,n)R n - 1 



In Table I are reported the expressions of a for some f(R) 
models. 

Notice that in the cases 1. and 2. of Table I, the 
quantity a is constant and independent of a, while the 
case 3. reduces to the case 2. as soon as 7 = 0. Moreover, 
the cases 3. and 4. yield a constant a in the regime of 
small or large scalar curvature. In general, a can be 
zero, negative or positive, depending on the values of the 
constants /3, 7, and n. 

The constants c(/3, n) and d(/?,n), appearing in Table 
I, are defined as 

c(/3, n) = 3(1 - 2/3) + (3/3 - l)n+ (111.5) 

+ M^i) [2/3 + 1 _ 2n]) 



d(/3, n) = 2/3 - 1 + n(l - /3) - 2n(n - 1) . (III.6) 

In what follows we shall focalize in particular on f(R) 
model given by the case 2. in Table I. 



IV. THE ABUNDANCE OF THERMAL RELICS 

According to standard cosmology and particle physics, 
the calculation of the relic density of particles is based 
on the assumption that the period of the Universe dom- 
inated by radiation began before the main production of 
relics and that the entropy of matter is conserved during 
this epoch and the successive one. Clearly, a different 
relic density of particles is expected once these assump- 
tions are modified. In this scenario, therefore, any con- 
tribution to the energy density modifies the Hubble ex- 
pansion rate, which reflects in a modification of the relic 
density values. Along these lines have been performed in- 
vestigations considering Brans-Dicke cosmological mod- 
els or anisotropic expansions (lfj, [ll| . 



The enhancement function A(T ) ap pearing in Eq. p.3[) 
is conveniently parameterized as [l(| (see also [3~lj ]) 

A(T) = { 1+71 (^j) T>T re (IV .l) 

I 1 T < T re 

Here T re denotes the temperature at which the Hubble 
rate reenters the standard rate of GR (to avoid con- 
tradictions with big bang nucleosynthesis, it is required 
T re > 1 MeV), while Tf is the temperature at which 
the WIMPs dark matter freezes out in the standard cos- 
mology (Tf ~ 17.3GeV). Notice that the value of Tf, 
in general, varies by varying the dark matter candidate 
mass m x . The parameters rj and v are free parameters 
characterizing the specific cosmological model. Estima- 
tions carried out in [10( have been obtained by setting 
T re = 1 MeV. If 77 ^> 1, the previous equation reads 

A(T)~Jpj U . (IV.2) 

The values of the parameter 77, required to explain the 
PAMELA data, are 

1 ^ V ^ !0 3 • (IV.3) 

The corresponding values of the WIMPs masses are in 
the range 

10 2 GeV < m x < 10 3 GeV . (IV.4) 

However, for dark matter masses of the order m x ~ 
10 2 GeV, the parameter -q can also assume values close 
to zero. Hence 

m x ~ 10 2 GeV < 77 < 1 . (IV.5) 

These considerations can be framed in the context of 
f(R) gravity where curvature energy density and pres- 
sure play a specific role in the characterization of the 
function A(T). 

V. THE AMPLIFICATION FACTOR A(T) IN 
f(R) GRAVITY 

We have now all the ingredients to study the role of 
modified cosmology provided by f(R) gravity to explain 
the PAMELA data via dark matter relic abundance. 

The — component of field equations (111.31) allows to 
write the relation between the expansion rates of the Uni- 
verse in f(R) cosmology and in the standard cosmology 
as 

where Hqr = (k 2 p m /3) 1 / 2 . Comparing with Eq. (PI) , 
one infers the amplification factor A(T) 

A(T) = Vm, (V.2) 
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where 



K 2 p m 



(V.3) 



The amplification factor A(T) strictly depends on the 
ratio between curvature and matter densities in the field 
Eq. (jV.f p and, as before discussed, we will deal with 
the regime where p c < p m , so that A(T) ~ 1 + | (in 
such a case A(T) is larger than 1 allowing to address the 
PAMELA result). 

The densities are, in principle, functions of the tem- 
perature T. The conservation of entropy S — constant 
implies that the scale factor and the temperature of the 
Universe are related by a = T /T [H (T = 3K is the 
present temperature of photons and we set oo = 1, where 
ao is the present value of the scale factor). The combi- 
nation of Eqs. (|II.16[) and (|II.14[) allows to rewrite Eq. 
(|V.3[) in the form 



K 2 P ™ 



P% (Tf 
Using the fact that 

rri 4 



3((7 — w) 



T,) 



3(ct— w) 



P a = pa = 



* 9* rj,i 



30 



J- n 



where are the relativistic degrees of freedom, and that 



ho^Orad — 6.15 X 10 , ^Orad — ^0 

finds 



0v 



Po = VoradPcr = 5.16 x 10- 51 GeV 4 , 

where p cr = 8.4/iq x 10 _47 GeV 4 is the critical density. 
Since Pq must be compared with K 2 p™, it is worthwhile 
to report the order of the magnitude of K 2 p™: 

« 2 p™ = 1.3 x 10~ 87 GeV 2 . 



In the regime r < 1 (A ~ 1 + r/2), Eq. (jlV.ip yields 
the following expression for the parameter v and r\ 



v = 3(cr — w) , rj 



Po 



2k 2 pS 1 \T q 



T 



3(cr — iu) 



(V.5) 



In the standard cosmological model, the Einstein field 
equations admit, for a radiation dominated Universe, the 
solution Hgr = l/2t, i.e. a = a^t 1 / 2 . On the other hand, 
if the time evolution of the cosmological background is 
governed by field equations (|II.3|) . then the scale factor 
a r~i t& yields 



H = t = 2pH GR -> p= + 



(V.6) 



Since (3 is constant (and positive), the r.h.s. of (|V.1[) 
(and Eq. (|V.4|| ) must be necessarily independent of the 



cosmic time, and therefore on the temperature. Such a 
condition is fulfilled by requiring 1 

v = <4 a-w = 0. (V.7) 

In literature, the case v = refers to an overall boost of 
the Hubble expansion rate. 

For r < 1, the constant f3 turns out to be related to 
the parameter r\ by the relation 



1 



+ V 3 
— <4 



as n — 



Po 1 
2k 2 p™ 2 



(V.8) 



However, constraints provided by the big bang nucleosyn- 
thesis give a more stringent bound. In fact, in the frame- 
work of f(R) cosmology one gets [36[ 



10" 



(V.4) so hat 



t) < 2 x 10" 



(V.9) 



(V.10) 



See also [32|, |33j for a general discussion on big bang 
nucleosynthesis . 

The previous results have a general validity because 
they do not refer to a specific form of f(R). As a conse- 
quence, f(R) cosmology represents a suitable framework 
for getting the amplification of the expansion rate of the 
Universe, and therefore to explain the PAMELA data. 
From Eq. (|IV.5I) , we therefore conclude that the mass of 
the relic dark particles is ~ 10 2 GeV. 

In what follows, we shall apply the above results to the 
model given by 



f(R) =ujR + aR n , 



(V.ll) 



where we insert the constant u to figure out terms coming 
from GR (at the end it will be set equal to 1). 



1 A comment is in order on Eq. tV.7t . Since our results are not 
exacts, we should more properly require v 0, i.e. a — w fa 0, 
with w = 1, a = w — g, and q <C 1 This means a very tiny 
variation of r with T. Another possibility to get r constant is 
to set a = w = ^ — Interactions among massless particles, in 
fact, lead to running coupling constants, and, hence, the trace 
anomaly T m oc /3(g)F^ v ^ 0. In Ref. [H it has been studied 
the thermodynamical potential of a plasma for SU(Nq) gauge 
theory, with coupling g and Nf flavors. This study shows that 
the adiabatic index (for the radiation) is given by w = i — q, 
where 



C 



f) 



18tt 2 (4tt) 2 



2+ 7_N c N f 



2Nl 



1 



up to 0(g 5 ) corrections (it is worth noting that typical gauge 
groups and matter content at very high energies can yield 1 — 
3w ~ 10" 2 - 10" 1 [3). 
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Let us start by checking the consistency between Eqs. 
(|lLT0l) and pill)) for this model. From Eq. ([HTTP]) , it 
follows 



(a,-l) 



(3R 



2(2/3-1) 



(V.12) 



aR 71 



2(2/3-1) 



[(2-n)/3-(2n-l)(n-l) 



a 2n//3 ■ 
Here p% is defined as 



(V.13) 



with 

r = 



2(2/3-1) 7J 

(V.14) 

The last expression has been obtained for u = 1. In 
the above equations we used i? = —6/3(2/3 — l)/t 2 and 
* = (a/ao) 1 ^- 

By comparing (|V.13j) with 15j) . one gets 



13 = 



2n 



3(<7+l) 



(V.15) 



c(/3, n) 

It is easy to verify that using a — — — - — - (case 2. in 

3a(/3,n) 

Table I), with c(/3, n) and d(/3,n) defined in (|III.5|) and 

(|III.6p , respectively, Eq. (|V.15|) is automatically fulfilled. 

In the case in which a « w « 1/3, Eq. (|V.15|) gives 

/8«^. (V.16) 
Equations (jESJ), (|VlO| and (|Vl6l) imply 



™ = 1 + »7 < 1 + 2 x 1CT 3 , /3~-. (V.17) 



and r\ assume the form 

2 ' ?/ U, 

(VI.18) 

The constant /3 turns out to be related to the parameter 
r\ by the relation 



2 ~ 2 



which means that the relic dark particles have masses 

> 10 2 GeV, see Eq. (jE SJ). 

Combining (|VL19|) with (|VTl"6|) it follows that in the 
regime r > 1 , the parameter n has to satisfy the condition 



n = Vv > 1 



(VI.20) 



Therefore, the model (|V.llj) solves the PAMELA puzzle 
as soon as the exponent n assumes values greater than 1, 
while the Universe evolves as a ~ t , with /3 > 1/2. 

Clearly, looking at the Table I, in the regime R 3> 0, 
the parameter a is a constant whose signature depends 
on the values of (3 and n. For example, considering the 
Case 3. in Table I gives 



/3 2 



1 



13 \P + 10 



(VI.21) 



This means that, for large /3, a approaches to -1 and then 
the curvature fluid behaves as a cosmological constant 
(see Eq. (|II.15[l ). In conclusion, it is possible to show that 
relaxing the BBN-constrain, i.e. during the pre-BBN era, 
there are several modified gravity models (in particular 
/(i?)-models) capable of fitting the PAMELA data. 



VII. EFFECTIVE POTENTIALS IN THE 
EINSTEIN FRAME 



Hence, the model works provided the constraint (|V. 17|) 
is satisfied. Such a result agrees with Solar System con- 
straints [35[ and the range in Eq. pV.5l) . 



VI. ANALYSIS IN THE PRE-BBN EPOCH 

Let us now study the case in which the BBN-constraint 
(|V.9[) is relaxed. This means that we are well before the 
BBN-cpoch. This case isiustified by the fact that, in 
the model proposed in [10(, the enhancement function 
A{T) is set to 1 as soon as the Universe reaches the BBN 
era. From this instant, it is assumed to evolve according 
to GR. In this perspective, it turns out interesting to 
analyze the case r > 1 where cosmological models are 
less fine-tuned. 

We can take into account the case r > 1 where A ~ 
y/r/2. From Eq. (|IV.2|) one gets that the parameters v 



The further gravitational degrees of freedom coming 
from f(R) gravity can be figure out as additional scalar 
fields. In this case, the enhancement necessary to explain 
the relic abundance (i.e. the amplification factor A(T)) 
can be tracked by the scalar field dynamics. 

Let us discuss the conformal transformation of f(R) 
gravity which can be transformed into Einstein's grav- 
ity plus a minimally coupled scalar field [H @ . We can 
assume that the energy density p c and the pressure p c 
are dominant with respect to matter. This assumption is 
very natural since curvature components can play the role 
of cosmological dark energy which is larger than matter 
component [H}. Furthermore, the assumption is justified 
by the fact that we are in the regime r > I. To better 
appreciate the role of the scalar field when passing to 
the Einstein frame, we can neglect its coupling to matter 
fields since it is not particularly relevant in this epoch 
and its signature can be hardly observed. 
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Let us consider the conformal transformation on the 
metric g^ v 



By choosing 



and setting 



— 6 ^9fJ>u ■ 



X =-ln|/'(i?)| 



kip = X ■ 



V6 



(VII.22) 



(VII.23) 



(VH.24) 



it can be shown that the Lagrangian density of f(R) in 
(11.11) can be recast in the (conformally) equivalent form 



-gf(R) 



while the field equations reads 



(VII.25) 



G>„ — K 



1 



V M V„ - -^g^vXI p ipXI p ip + g^V(ip) 



The potential V is defined as 

f-Rf 



v = 



2f 



(VII.26) 



(VII.27) 



Let us figure out the explicit form of the potential V in 
the case in which f = R + aR n . Inverting (|VII.23[) , one 
gets 



/ _ 2ktp 



(VII.28) 



Using /' = 1 + anR n 1 , the previous equation yields 
1 



R 



-{e 2kip - 1) 



The potential (|VII.27|) reads 



V 



a(l - n) ^_ 4kip 



an 



(VII.29) 



(VII.30) 



2-'tt(l-") c it^ 
(cm) "- 1 



[sinh kip] 



For kip <^ 1 the potential assumes a power law behavior 



where 



V ~ V tp d , 



_ 2"-i q(l-7i) _ 7i 
VO = : — , = 



(VII.31) 



(an) ' 



n-1 



Such a form of potential has been widely studied in liter- 
ature in the framework of alternative theories of gravity 
1- 




k (p 



FIG. 1: V/Vo vs kip for n = 2. V/V -> 1/4 as k<p > 1. 




k i£ 



FIG. 2: F/Fo vs kp for n = -0.5, n = 1.5, and n = 3. for 
large kip, the potential approaches to zero for n = —0.5 and 
n = 3, while grows up for n = 1.5. 



It is worthwhile to point out that, in the case n = 2, 
the potential (|VII.30P tends to a constant value for large 

kip 



V 



4 



as kip^> 1 



(VII.32) 



Therefore, in this regime the potential plays the role of 
a cosmological constant. This behavior is shown in Fig. 
[TJ where it is plotted V/Vq vs kip for n — 2. To compare 
with other values of n, we also plot V/Vq for n = —0.5, 
n = 3 and rt = 1.5, see Fig. [2j In the first two cases, the 
potential V approaches to zero as kip 3> 1, in the latter 
V grows for increasing kip. 

Specifically, the amplification factor (|V.2[) can be ex- 
pressed in terms of the scalar field ip. By using the scale 
factor (|III.1|) and the relation between the Ricci curva- 
ture and the scalar field, Eq. (|VII.29|) . one obtains 



afeO = ao|6)8(20-l)|*(an) 



- : - 2l^T) ( e 2fe¥> 



(VII.33) 

As a consequence, the amplification factor A(ip) reads 

A(ip) = [l + r(ip)] 1 / 2 , (VII.34) 
where the quantity r, defined in (|V.4|) . is given by 



iip)=r a [2 2 ^-l]" 



(VII.35) 
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with 



ro 



Po 



00160(20-1 







(VII.36) 

These considerations mean that it is straightforward to 
translate results in Sec.V in the Einstein frame and the 
amplification factor, being linked to the Hubble param- 
eter, strictly depends on the form of the potential V((p) 
and then (p. In other words, dynamics related to p c and 
p° can be recast in terms of V(tp) and ip. 

Reversing the argument, the PAMELA data (an then 
the requested amplification factor to explain them) could 
be a probe for alternative theories of gravity, in particular 
for f(R) gravity. 

Finally, as discussed in (32H29J, it is worth noticing 
that f(R) gravity can, in general, give rise to cosmologi- 
cal viable models compatible with radiation and matter- 
dominated epochs evolving into late accelerated phases. 
This means that reliable models can be phenomenologi- 
cally reconstructed by means of observational data. 



The analysis performed in this paper relies on f(R) 
models assumed as function of the scalar curvature R. 
However, other curvature invariants like Riemann and 
Ricci invariants can be considered to refine the analysis. 
As final comment, it is important to stress that PAMELA 
data are very strict results indicating that revisions are 
necessary at fundamental physics and cosmological levels. 



Appendix A: Power law solutions, continuity and 
field equations 

In this Appendix we provide details of formulas used 
in the main text. For the FRW metric and for a scale 
factor evolving as a power law, a = aot" , we have 



H 



Rq 



-3(H + H 2 



o _ 0(0 - 1) 
*' a P ' 

30(0-1) _ 1 0-1 
i 2 ~ 2 20 - 1 



(A.l) 



R, (A.2) 



VIII. CONCLUSIONS 

There is nowadays a great interest toward the recent 
data of PAMELA experiment since they could represent 
a possible signature for dark matter. Among the vari- 
ous mechanisms that have been proposed to explain the 
PAMELA results, the models based on alternative cos- 
mologies represent suitable candidates. These models 
rely on the fact that allow an enhancement of the early 
expansion rate of the Universe, hence large annihilation 
cross sections. At the same time, such models are also 
compatible with other observations [J- 

In this paper, we faced the problem of the PAMELA 
puzzle in the context of cosmological models provided 
by f(R) gravity. Our approach is based on the effec- 
tive description of the energy-momentum tensor (the to- 
tal energy-momentum tensor is written as the sum of 
the energy- momentum tensor of ordinary matter, T™, 
plus the effective energy-momentum tensor T° induced 
by curvature terms appearing in the nonlinear action 
of gravity - these two fluids satisfy separately Bianchi's 
identities) and on the assumption that the scale factor 
evolves as a power law (a ~ t@). By confining ourselves 
to the case in which T™ > T° , we found that value of the 
dark matter masses necessary to explain the PAMELA 
puzzle turn out to be of the order m x ~ 10 2 GeV. In 
particular, for the model f(R) = R + aR n it follows 
that to explain the PAMELA data, the admissible value 
of the exponent n is given by Eq. (1V.17[) at the BBN 
epoch. Relaxing such a hypothesis and considering also 
the pre-BBN epoch, the model is less "fine-tuned" and 
other f(R) models agrees with data. The evolution can 
be tracked considering conformal transformations. In the 
Einstein frame, the amplification factor A depends on a 
scalar whose potential determines the model dynamics. 



R=-«H + 2H* ) = - W »- 1) , (A.4) 
R=-&{H + AHH) = -^-, (A.5) 

R = -6(H + 4iJ 2 + AHH) = ^ , (A.6) 

□/' = /' + 3Hf' , (A.7) 

VoVo/'-/', VX l f ///'-v;. (A.8) 

Let us now show that Eq. (|IL1~2"]) is fulfilled assuming 
that the Universe evolves with a scale factor of the form 
a ~ t p . By using the relation R = -2R/t = -2HR//3, 
one can write ff as 

p c = -^^(-l + f-f"R)-3Hf', (A.9) 

On the other hand, summing up (|III.2[) and 3() and 
using HR = i? 2 /3(20 - 1) it follows 

P c +P C = 3(2/ f_ 1} (-1 + /' - f"R) + f ■ (A.10) 
These equations show that Eq. (|II.12|) is fulfilled. 



9 



Finally, we discuss in some details the solutions of the 
field equations (jll.ip in the regime r < 1. For the (|II.9|) 
they read 

-3-f-£ + 3-f"R = « l tf n , (A.ff) 
a 2 a 



3/3 



6(2/3 -f) 



ujR 



ac(f3,n) 
6(2/3 -f) 



iT = K 2 W P V 



(A.14) 



a a 
- +2 

a 



/ 



, , f + J - - 2-f'R - f'R 2 f'R = «V 
/ 2 a 



(A.12) 

For the model (IV.llj) and a = a t^ , the above equations 
can be recast in the following form 



2/3 



2/3- 1 



loR 



ad{P,n) _ 2 
2(2/3 -F) ~ KP 



(A.13) 



where we have used the above results and p m — wp m 
(w = 1/3), while c and d are defined by Eqs. (|III.5[) and 

CBDp . 



In the regime r < 1, therefore the i? ra -term is a per- 
turbation, Eqs. (IA.13|) and (IA.14|) show that a possible 
solution is provided by n ~ 1 , according to (|V.17I) . 
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